Gravitational waves from a spinning particle scattered by a relativistic star: 

Axial mode case 
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We study gravitational waves from a spinning test particle scattered by a relativistic star using a 
perturbation method. The present analysis is restricted to axial modes. By calculating the energy 
spectrum, the waveforms and the total energy and angular momentum of gravitational waves, we 
analyze the dependence of the emitted gravitational waves on a particle spin. For a normal neutron 
star, the energy spectrum has one broad peak whose characteristic frequency corresponds to the 
angular velocity at the turning point (a periastron). Since the turning point is determined by the 
orbital parameter, there exists the dependence of the gravitational wave on a particle spin. We find 
that the total energy of I = 2 gravitational waves gets larger as the spin increases in the anti-parallel 
direction to the orbital angular momentum. For an ultracompact star, in addition to such an orbital 
contribution, we find the quasi-normal modes exited by a scattered particle, whose excitation rate 
to gravitational waves depends on the particle spin. We also discuss the ratio of the total angular 
momentum to the total energy of gravitational waves and explain its spin dependence. 

PACS number(s): 04.25.Nx, 04.30.-w, 04.40. Dg 



I. INTRODUCTION 

A neutron star binary is one of the most promising 
source for gravitational waves, which will be observed 
in the near future by gravitational wave interferometers 
such as Laser Interferometric Gravitational Wave Obser- 
vatory (LIGO), VIRGO, GEO600, TAMA300, and Laser 
Interferometer Space Antenna (LISA) |l|. Since gravi- 
tational waves has a little interaction to the matter, we 
could resolve the phenomena of a strong gravitational 
field such as a black hole formation by the direct obser- 
vation of gravitational waves, which can not be revealed 
completely only by electromagnetic waves. From such an 
analysis, we could find new information about the inte- 
rior of a neutron star whose equation of state is not yet 
well-known. These expectations motivate theoretical ef- 
forts to investigate gravitational waves from a neutron 
star binary in the inspiral and merging phase. 

One of the final step to calculate gravitational waves 
from a neutron star binary is a fully general relativistic 
simulation with no approximation. In the past decade, 
many efforts have been put into the 3D numerical sim- 
ulation of merging a neutron star binary in Newtonian, 
post-Newtonian, and general relativity. However, Ein- 
stein equations are non-linear wave equations with con- 
straints. The present simulation has still faced to many 
difficulties, however we hope that we will succeed to sim- 
ulate some of limited cases in the near future. In such a 
situation, we may invoke another way to study gravita- 
tional waves from such a binary system by using some ap- 



proximation methods focusing to some particular phys- 
ical quantities and extract fundamental and important 
properties. 

One of these approximate approaches of gravitational 
waves from a binary system is a perturbation method, in 
which a companion of the binary is treated as a test parti- 
cle with mass /z moving around a relativistic background 
object with mass M (^ /i). It has been shown that this 
approach gives a good approximation for gravitational 
waves from a head-on collision of two black holes [EJ , and 
it has been thought to be one of the reliable approxi- 
mation for treating gravitational waves from a coalescing 
black hole binary. 

This method can also apply to a neutron star binary, 
although so far there has been a strong restriction such 
that we assume a spherical symmetric background space- 
time, which can be extended to a slow rotating back- 
ground. The motion of a test particle excites nonradial 
pulsations of a central spherically symmetric star in the 
emission of gravitational waves. The perturbation equa- 
tions of gravitational waves outside the spherical star are 
given by Regge and Wheeler |0] and Zerilli H|. Regge- 
Wheeler-Zerilli equations are well-known as perturbation 
equations in the Schwarzschild background spacetime. 
While inside the star, the perturbation equations both 
for the metric and stellar fluid were derived by Thorne 
and Campolattaro |q], Lindblom and Detweiler IgI, and 
Chandrasekhar and Ferrari . 

For a spherically symmetric star, the perturbation 
equations are completely separated by its parity, i.e., po- 



lar (even) and axial (odd) modes. The polar mode con- 
tains two types of waves; gravitational waves and mate- 
rial waves which interact each other. While, the axial 
mode has only one type, i.e., gravitational waves. From 
this character, quasi-normal modes of fluid oscillations 
(/, p and g modes) are exited by a test particle only in 
polar modes, while the quasi- normal modes of gravita- 
tional waves {w modes) will appear in both modes. As 
for the axial mode, since there is no fluid oscillation, only 
w modes are excited. 

In the case of a binary black hole system, it is known 
that gravitational waves are dominated by the quasi- 
normal modes at the ringing-down phase. Since the 
quasi-normal mode contains many information about the 
newly formed black hole or neutron star, the analysis of 
the quasi-normal modes is important in the relativistic 
astrophysics. As for a binary neutron star system, the 
remnant of merged binary stars with M ~ IAAIq may 
form a black hole after a short time. The ringing-down 
phase of a neutron star binary may be determined by the 
quasi-normal modes of a black hole. However, since the 
quasi-normal modes of a star would also be excited at 
the inspiral phase, it is still important to focus on them. 

For polar modes, Kojima |q| examined gravitational 
waves from a particle in a circular orbit around a spher- 
ically symmetric relativistic star. From the symmetry of 
a circular orbit, axial modes are completely canceled out 
in this case. He found the resonant oscillation of a star 
and an enhancement of gravitational waves at the same 
frequency of the quasi-normal mode. 

For axial modes, Borrelli |^ calculated the energy spec- 
trum when a particle is spiraling onto an ultracompact 
star {R < 3M). In the case of the ultracompact star, it 
is known that some quasi-normal modes are trapped by 
a potential barrier of gravitational waves. These quasi- 
normal modes are excited by a falling particle. 

Several works for a scattered particle have been done. 
The present authors |lO[ calculated the energy spectrum 
for axial modes. In the case of neutron stars {R > 3A/), 
a single peak appears in the energy spectrum, which de- 
pends on the trajectory of a particle. On the other hand, 
in the case of ultracompact stars (i? ^ 3M), many sharp 
peaks are also excited in the energy spectrum. Some 
peaks are explained by the quasi-normal modes trapped 
by the quasi-bound state in the potential, and the rest of 
peaks is due to the resonance of two propagating waves, 
which are reflected by the potential barrier or the infinite 
wall at the center of the star. These peaks never appear 
for the black hole case, because there are no bound state 
in the potential, and gravitational waves are purely in- 
coming at the horizon, so the resonance of two reflected 
waves never occurs. Andrade and Price |ll|| also stud- 
ied axial modes of gravitational waves by a scattered 
particle. They have found that even for stellar models 
with R > 3M, a scattered particle with a very large 
orbital angular momentum could weakly excite quasi- 
normal modes. Ferrari, Gualtieri, and Borrelli ^^ com- 
puted energy spectra of both polar and axial modes from 



a particle scattered by a polytropic star. In their results, 
fluid modes (/ and p modes) are also excited in the en- 
ergy spectrum. 

Recently, there are two papers demonstrated gravita- 
tional waves from a scattered particle in a spherical sym- 
metric star for time evolution. Ruoff, Laguna, and PuUin 
p3[ computed spectrum and waveform for polar modes, 
focusing on the excitation of lu-modes. They concluded 
that the excitation appeared in the spectrum, although 
w-modes significantly contributes to the energy in the 
case of the orbital speed v > 0.9c. Ferrari and Kokko- 
tas p4| also computed spectrum and waveform for axial 
modes. They also point out that most of the w-modes 
are indeed, clearly excited in the spectrum whose feature 
is the same to our previous results |lO|| . 

Neutron stars are usually rotating. Then it is very 
natural to take into account a rotation effect as a next 
step. From the viewpoint of the gravitational wave as- 
tronomy, a rotating relativistic star is a very important 
object, because there exists r mode which shows very 
strong instability for any rotations ||l^ and this effect it- 
self becomes one of the promising source of gravitational 
waves pq ]. For example, hot young neutron stars would 
be enough rapidly rotating to occur r mode instability, 
and would spin down toward a usual observed neutron 
star by radiating gravitational waves. The recent obser- 
vation suggests that low mass X-ray binaries p7| would 
include a rapidly rotating neutron star. 

Although there are several perturbation analysis for 
a slowly rotating neutron star, we have no systematic 
method to solve for the case of a rapidly rotating star 
except for the analysis of the neutral point |1^. Here, we 
propose an extremely simple model to treat some effect 
of a rapid rotation by adding a spin to a test particle. 
Since we can not take into account the structure of a test 
particle, we do not expect the excitation of a test particle 
itself, and then will not be able to discuss the wave mode 
of a spinning particle itself. However, this work might 
be a milestone for a further study of a rotation effect on 
gravitational waves from a binary system. 

The spinning particle in a curved spacetime was for- 
mulated by Papapetrou ||l^ and Dixon |20J by taking a 
point particle limit of a relativistic extended body. In the 
case of a rotating black hole with a spinning test parti- 
cle, which mimics a collision of two rotating black holes, 
there are several works using a perturbation approach. 

In this paper, we will examine gravitational waves from 
a spinning particle scattered by a relativistic star. We 
calculate the energy spectrum, the waveform and the to- 
tal energy and angular momentum of the emitted grav- 
itational waves, and investigate the spin effect to grav- 
itational waves. This paper is organized as follows. In 
Sec. ||, we briefly review the equation of motion of a 
spinning particle and gravitational waves in a spherically 
symmetric background using a perturbation theory. We 
show our numerical results in Sec. 
devoted to discussion. 



[II. Section [V is 



Throughout this paper, we adopt the units of G = c = 
1, where G and c denote the gravitational constant and 
speed of hght respectively, and the metric signature of 



II. PERTURBATION EQUATIONS FOR AXIAL 
MODES 

A. Spinning Particle 

We briefly summarize the equations of motion of a 
spinning test particle and its energy-momentum tensor 
T^'^ . The motion of a spinning test particle in a curved 
background is discussed by Papapetrou ||I^ and Dixon 
pO| . They presented the equation of motions for a spin- 
ning particle as 



^P^{t) = -\r^^.,Az{tW{t)SP^{t), 



D 



5'^'^(r) ^p^'{T)v''iT)~p''{T)v^'{T), 



(2.1) 
(2.2) 



where z''{t), v''{t) = dzf^/dr, p^'{t) and S'"'{t) denote 
the trajectory of a particle, 4-velocity, 4-momentum, and 
spin tensor, respectively. The mass ^ and the magni- 
tude of a spin S, which are defined by p^^p^ — — /i^ and 
SfivS^^ = 25*^, are conserved [Q. Since we also have to 
specify the center of mass, which supplies an additional 
condition, we assume 



5^>. = 0, 



(2.3) 



with which we find a complete set of equations of mo- 
tion. The normalization of the affine parameter r is con- 
veniently chosen as m^w'^ = —1, where u^{= p^ j [i) is the 
normalized momentum. 

Since the total angular momentum of the particle is 
conserved, we can specify the z axis (and then the equa- 
torial plane) in Schwarzschild spacetime by the direc- 
tion of the conserved total angular momentum J, i.e., 
J — (0, 0, J) with J > 0. In this paper, we restrict 
a particle motion onto the equatorial plane (Q — 7r/2) 
in Schwarzschild background. By use of the conserved 
quantities pfl] along the trajectory z^(t\ we write down 
the equation of motion of the particle with the energy 
i5(= E / ^), total angular momentum J(= J/ ^) and spin 
S{= S/fi) H as follows ||: 
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dt 
d^ 



2M 



,^E-Mli\, p.4) 
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(2.5) 
(2.6) 



S.s = r^ 1 



MS^ 



3MS^r 



. _\ 2 



Ac = 1 - 



( J - SE^ 



y3 



(2.7) 
(2.8) 
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^ ^-(--) 



(2.9) 



In what follows, we introduce the "orbital" angular mo- 
mentum by L = J— S" which is same to Ref. |2J , because 
the orbital angular momentum plays an important role 
in the gravitational wave in the case of a non-spinning 
particle [Q. Since we set J > 0, the spin direction is 
classified by the sign of S into the parallel {S > 0) and 
anti-parallel (5 < 0) to the orbital angular momentum. 
From Eq. (pj|), we can introduce the effective po- 
tential y(particio) Qj: ^ spinning particle in the equatorial 
plane, which is 



^(particle) ^ -/? + V (3^ " ^7 

a 
where a, /3, 7 are 



r'-S'll- 



2M 
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(3:^r{r- 3M)SJ 
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-T) 



- 1- 



2M 



(2.10) 

(2.11) 
(2.12) 

(2.13) 



where S5, A5 and Rs are 



Since we have not normalized the 4-velocity v^ , we may 
find some unphysical orbital motion with v^'^ > for 
some specific choice of orbital parameters p3[ . In order 
not to take such unphysical parameters, we impose the 
timelike condition, v^^v^ < 0, which is shown in Fig. n. 
The dotted line denotes a null condition, i.e. Vf^v^ = 0, 
and the timelike condition is satisfied in the region above 
the dotted line. The solid line denotes the maximum 
point of the effective potential y (particle) jj- ^ particle 
has smaller energy than that on the solid line with a 
given L, it will be scattered by a potential barrier, which 
we are studying here. The forbidden parameter range 
appears as the particle spin gets larger. 

In order to study the emitted gravitational waves from 
a test particle around a relativistic object, it is important 
to specify the turning point, i.e. a periastron. Because 
the turning point corresponds to the strongest gravita- 
tional field in the particle orbit and its angular velocity 
gives a characteristic frequency of the emitted gravita- 
tional waves 1 1^1 . In Fig. 0, the location of the turning 
point rniin and its angular velocity d(j)/dt are shown in the 



case of the normal neutron star with R = 5.QM. We find 
that Tmin increases (and then dip/dt decreases) monoton- 
icaUy as the spin S gets larger in the parallel direction. 
However, if the stellar model is an ultracompact star, e.g. 
with R — 2.26M, and a particle passes through near the 
ultracompact star, the feature at the turning point dras- 
tically changes, and there exists a minimal point for the 
location and angular velocity as shown in Fig. ^. This is 
just because the particle passes through a strong gravita- 
tional region, and then nonlinear coupling term between 
L and S will appear. 

The energy momentum tensor of a spinning particle is 
given ^ by 



T"" = / dT 



fu/ ^ i^)r sS'-^Kx- z{t)) 



\/p{SP^^'{x,T)v^\x,T) 



(5(4)(x~z(r)) 



-9 



, (2.14) 



where v'^{x, r), p'^{x, r) and S^'\x, r) are extended from 
the 4-velocity, 4-momentum and spin tensor by the par- 
allel displacement bi- vector gf^aix,x') [^ as 

v^ix,T)^g\{x,ziT))v''{T), (2.15) 

p^^{x,T)^r<.{x,ziT))p''iT), (2.16) 

S'^-'ix, t) = ro.{x, z(T))rp(x, z{t))S'^P{t), (2.17) 
where the bi-vector ^''^(x, z(t)) satisfies 

lim5^"(x,z(r))=<5^", (2.18) 

limV,g^„(a;,z(T)) = 0. (2.19) 



B. Linearized Einstein Equation 

We discuss gravitational waves from a spinning test 
particle scattered by a relativistic star. We assume that a 
star is spherically symmetric. The spherically symmetric 
background metric is described as 

g^^}dx''dx'' = -e''('')dt2 + e^^'^'dr^ + r^ (^^2 _^ ^^^^2 ^^^2^ 

(2.20) 

We shall introduce a mass function M{r) as 

2M{r) 



-Hr) 



1 



(2.21) 



The gravitational mass of a star is given by M = M{R) 
where R is the surface radius. We also assume that the 
background star is composed of a perfect fluid 



Tf"" = {p + P)u''u'' + Pg'"', 



(2.22) 



where p and P are the density and the pressure of the 
fluid, respectively. 



In a spherically symmetric background, we can al- 
ways decompose the perturbed metric into an axial mode 
l^(axia ) ^^j ^ polar mode hfF"'^^ . In what follows, we 
discuss only an axial mode. In the previous paper [[lO| , 
we presented the formalism to calculate the gravitational 
wave from a test particle scattered by a star. In this pa- 
per, we consider a spinning test particle and discuss the 
effect of a spin. Since the formalism is the same except 
for the source term, which depends on the character of 
the particle, we present a brief sketch of the formalism. 
We use the same notation and variables as those in Ref. 

It is known that, using the Regge- Wheeler gauge, per- 
turbation equations of axial modes are reduced to a single 
wave equation both in the interior and exterior regions of 
a star. The perturbation equation inside a star is given 

by 



d^x. 



(int) 



dr*2 



(c^ - vr') 



X. 



(int) 



= 0, 



(2.23) 



where the "tortoise" coordinate r* and the effective po- 



tential Vi ' are defined as 



-(-A)/2^ 



V, 



(int) 



l{l + 1) 6A/(r) 



47r(P - p) 



(2.24) 



(2.25) 



respectively. 

The perturbation equation outside a star is derived as 



,2 -i^(cxt) 

dr*^ 



+ (lo^ - y/^"*^) X 



(cxt)^ Tj^(cxt) _ (,(cxt) /„ „„N 



where the tortoise coordinate r* and the Regge- Wheeler 
potential V^ are given by 



V' 



.*=.-, 2Mln(^-l), 
(e.t) ^L_ 2M\ f l{l + \) _ 6M 



r / \ r 



(2.27) 
(2.28) 



respectively. The source term Sj^^J in Eq. ( 2.26| ) is given 



by the ene rgy-m omentum tensor of a spinning test par- 
ticle (Eq. ( 2.14 )) expanded by tensor harmonics, whose 
explicit description is presented in the Appendix (see also 
Eq. (2.24) in Ref M). 



To construct the wave functions, X 



(int) 



and X 



(cxt) 



we have to impose a set of boundary conditions. The 
wave function Xi™ inside a star should be regular at 
the center of a star, which implies 



X 



(int) 



= m^'^ 



l+l 



1 



2(2/ + 3) 



47r(r 



+0{r^) 



2)(-(2/-l)pe-Pc 



— uj'^e 



(2.29) 



where rjii^ is an arbitrary constant, and Pc, Pc, and Vc 
are the central values of the density p, pressure P, and 
metric function v, respectively. For the wave function 
-^imw outside a star, there are no incoming waves at 
infinity, which is guaranteed by the condition 



X^ 



(cxt) 
IniLO 






(r* -^ CX3), 



(2.30) 



where Ai^^^ is the amplitude of an outgoing wave at in- 
finity. 

With those boundary conditions, we impose the match- 
ing condition at the surface r* = R* as 



x[T\R*)-xt^l\R*l 



dX, 



(int) 



dr 



^(^* 



dX. 



(oxt) 



dr* 



(R* 



(2.31) 
(2.32) 



which guarantees that the wave function is continuous 
and smooth. In this way, we construct the wave function 

Finally, we introduce important physical quantities of 
gravitational waves, i.e., the total energy, the total an- 
gular momentum, the energy spectrum, the angular mo- 
mentum spectrum, and the waveform as 



e: 



(axial) 
GW 



J, 



(axial) 
GW 



(IT \ (axial) 



(2.33) 

(2.34) 






327r 



lil + l)il-l)il + 2) 



A 



(axial) 
Imuj 



f dJawV"'"'"'^ m f dEcwV"'"'"'^ 



\ d(jj 



J (axial) I . , (axial) 






(2.35) 
(2.36) 



(axial) 
rauj 



l.rn 



X V;(^ + l)(^-l)(? + 2)e-""(*-''') 
x±2Yira{e,(l)), (2.37) 

where ±2^m is the spin-weighted spherical harmonics. 



III. GRAVITATIONAL WAVES FROM A 
SCATTERED SPINNING PARTICLE 

We analyze the emitted gravitational waves from a 
spinning test particle scattering around a relativistic star. 
Because the potential of the linearized Einstein equation 
takes a maximum value around r ~ 3Af , our stellar model 
can be classified into two models: One model is a nor- 
mal neutron star (i? > 3M) whose potential decreases 
monotonically as the coordinate radius increases. The 



observed neutron stars are classified into this case. The 
other model is an ultracompact star (i? < 3M). It would 
be very exotic, but show many interesting relativistic ef- 
fects, e.g. w modes trapped by a quasi-bound state of 
the potential (see, for example Figs. 5 and 6 in Ref. 
|lO|). We only discuss most simple equation of state, i.e., 
a uniform density star {p = const.), because the fluid 
oscillation will not couple directly to the axial modes of 
gravitational waves and then the information of the stel- 
lar matter may not be important. 



A. Normal neutron star 

Here, we discuss a normal neutron star with R ^ 3M. 
In the previous analysis of a non-spinning particle |1G| , we 
showed that one broad peak appears in the energy spec- 
trum of gravitational waves, which corresponds to the an- 
gular velocity of a test particle at the turning point, i.e., 
a periastron. Since the orbit of a test particle is affected 
by its spin, we expect that the existence of a spin will 
modify the spectrum of the emitted gravitational waves. 
In fact for the case of a spinning particle (Fig. ||), we 
find a broad spectrum quite similar to the non-spinning 
case, but we see that the energy spectrum is affected by 
a particle spin. The force which acts on the particle is 
given by 



/^(particle) 



M 



ni^(particle) 



dr 

M 
r 



1 



.(^^ 



NP 



") 



■{L^-ALS 



1 7- 



2M- 

(3.1) 

As we can see from this expression, the coupling between 
a spin 5" and an orbital angular momentum L, i.e. L ■ S 
coupling, gives an repulsive force. As a result, the turn- 
ing point leaves away from a stellar surface for a large 
positive (parallel) spin and then the angular velocity of 
a particle at the turning point gets smaller (Fig. 0). In 
fact, for a large positive (parallel) spin, a broad peak in 
the spectrum shifts to the lower frequency region, and for 
a large negative (anti-parallel) spin it shifts to the higher 
frequency region (Fig. ^). The change of the spectrum 
strength also corresponds to the shift of the turning point 
location. The amount of gravitational waves increases 
when the particle gets a large negative (anti-parallel) 
spin, because the gravitational field becomes strongest at 
the turning point and the most amount of gravitational 
waves radiates at this point. 

The total energy of gravitational waves is shown in 
Fig. H for each multipole mode I. As we expect, the 



I = 2 mode is dominant and the amount decreases as 
I gets larger. As for the spin dependence, we find that 
for small multipole modes {2 < I < A) the total energy 
monotonically increases as the spin gets large negative 
(anti-parallel). While for large multipole modes (/ > 5) 
such a simple relation disappears, but the energy reaches 
some maximum value for some negative value of S. We 
would expect that there is a complicated coupling in the 
higher order of Eq. (3.1) and its effect might directly 
appear to the total energy of gravitational waves. 

We show the ratio of the total angular momentum to 
the total energy, i.e. Jgw/^GWj in Fig. ^ (a). Re- 
markably, we find a simple relationship that Jqw/Eqw 
is proportional to the particle spin, and depends less on 
I. For the case of the scattered particle, the characteris- 
tic frequency of gravitational waves is described by the 
angular velocity at the turning point, i.e., 



(3.2) 



a;eh = mi- 



We roughly estimate the total energy and the total an- 
gular momentum of gravitational waves as 



iEo.t-'^Jdu.j:{^) 
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(axial) 
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(axial) 



f (IEqy^ 



(3.4) 
where AtJ ch is the width of the broad peak. Using Eqs. 



( p. 3D and (|3.4D, we find the relation 

/ J \ (axial) / 7 ^ \ —1 

/ •Jgw \ 



\ Egw J I 



dt 



(3.5) 



In fact, our numerical result has a good coincident with 
this rough estimation (Fig. || (b)). We then conclude 
that the spin dependence on (Jgw/-E'gw)/ comes from 
the angular velocity at the turning point. 

As for the waveform, a particle spin does not make a 
drastic change of the global feature of it (Fig. W. The dif- 
ference depending on the particle spin is only the ampli- 
tude of the burst wave. Since the burst wave is generated 
when the particle comes to the turning point, the ampli- 
tude is more enhanced when the particle can move in the 
stronger gravitational field. This is the reason why the 
amplitude is larger when the particle has an anti-parallel 
spin with the fixed orbital angular momentum. 

We conclude for the case of a normal neutron star that 
the effect of a spin could be understood just by the shift 
of the turning point, which is mainly caused hy L ■ S 
coupling. 



B. Ultracompact star 

If a star is sufficiently small as an ultracompact star 
{R ;$ 3Af), t he e ffectiv e pot ential of the gravitational 
waves (Eqs. (2.25) and (2.28)) has a minimum and the 



gravitational waves could be trapped and enhanced as a 
quasi-normal mode. This provides an additional charac- 
teristic structure in the spectrum of the emitted gravita- 
tional waves. In the previous works for a non-spinning 
particle |lQ,Ol , we found that a test particle scattered by 
an ultracompact star excites many quasi-normal modes, 
and there appear many sharp peaks in the energy spec- 
trum, besides one broad peak corresponding to the turn- 
ing point. We also found that such excited quasi- normal 
modes will leak later from the potential barrier and the 
resonance waves appears periodically in the waveform. 
Therefore in the case of an ultracompact star, there are 
two aspects in the spectrum: One is a broad peak which 
depends on the particle orbit and the other is many sharp 
peaks which depend on the structure of the background 
star. Here, with respect to the effect of a spin on the 
above properties, we shall discuss each in order. 

If the turning point is far from a stellar surface, we find 
the similar spectrum to the case of a normal neutron star 
except for many sharp peaks. The energy spectrum is 
shown in Fig. H, in which a particle with E = 1.01/i and 
L = A.b^iM is scattered at r^m - 5.61M {S = -O.S^iM) 
and r,„i,i ~ 7.31A/ (5 — O.S^M) by an ultracompact 
star with R = 2.26M. As the case of the normal neutron 
star, the broad peak shifts to the lower frequency region 
for a large positive (parallel) spin, while to the higher 
frequency region for a large negative (anti-parallel) spin. 

We also find many sharp peaks, which correspond to 
trapped quasi-normal modes induced by a scattered par- 
ticle. Since the quasi-normal modes depend only on a 
star, not on the property of a test particle, the frequency 
of those peaks completely corresponds for any values of a 
spin. The difference appears only in the strength of these 
peaks (see the small figure in Fig. g, which is enlarged 
on the frequency region of 1st quasi-normal mode). In 
the high frequency region, however, we also find some 
peaks in the spectrum, which show a little dependence 
on the particle spin. This is because those peaks do not 
really correspond to the quasi-normal mode, but appear 
as a result of an interference effect between incoming and 
outgoing gravitational waves, which depends on a spin, 
reflected by the potential barrier around r ~ 3A/ or the 
infinite wall at the center of a star (see Ref . |1^ ) . How- 
ever, since the contribution of such peaks to the emitted 
gravitational waves is very small, we may not be able to 
distinguish a spin effect in the gravitational waves. 

If a test particle approaches very close to the stellar 
surface, e.g. the turning point is r,„in ~ 3.15Af ~ 1.39-R 
(5 = 0) and r,-nin ^ 3.49Af ~ IMR {S = 0.8/iA/) for 
a particle with E — 2.38/i and L = 12.0/iA/, we find 
qualitative difference in the spectrum (Fig. ||). Here, we 
note that in this case a test particle with an anti-parallel 



spin {S < 0) is not allowed to give a scattered orbit, 
and any particles with such parameters are falling onto 
a star (More precisely, the lower limit of the spin value 
to allow a scattered orbit is 5 ~ — 0.012/iM). Thus we 
discuss only the case of a parallel spinning particle. For 
a parallel spin (5 > 0), the energy spectrum of gravita- 
tional waves is given in Fig. M. As the spin increases, a 
broad peak in the spectrum is monotonically shifted to 
the lower frequency region, while the angular velocity at 
the turning point has the minimum point at some spin 
S ~ 0.3^M (see Fig. |l^. To explain this broad peak 
frequency, the turning point is no longer a good indica- 
tor, although it is still one of the important factor to 
describe the emitted gravitational waves. Since the turn- 
ing point is very close to the surface of an ultracompact 
star, we suppose that there is another factor relating to 
the effect of a strong gravitational field, which probably 
affects to gravitational waves. In fact, the location of a 
turning point comes near to that of the maximum point 
of Regge- Wheeler potential in this case. 

Since quasi-normal modes are excited by a test parti- 
cle and trapped in the potential well, those modes will 
leak periodically after the encounter. In fact in Fig. |I|, 
we find that after the burst wave the ringing waves are 
emitted first a,t t — r* ^ lOOAf , and those are repeated 
with a smaller amplitude than the preceding one. Since 
quasi-normal modes are independent of the character of 
a particle, we do not find any dependence of the spin in 
those ringing phase except for its amplitude. 

The total energy of the emitted gravitational waves 
are shown in Fig. O for a particle with E = 2.30/x and 
L = 12.0/iM. Here we show the total energy only for 
I = 2,3, because if the multipole mode I becomes higher, 
some of the imaginary part of quasi-normal modes be- 
come extremely small. The existence of these small imag- 
inary parts requires to use very small step size of the 
frequency in order to integrate the energy spectrum pre- 
cisely, and then the analysis becomes quite difficult and 
time consuming. We also believe that those two multi- 
pole modes are dominant and then the total energy will 
not change so much (e.g. for the normal neutron star in 
Fig. |[ the energy of Z = 4 is about 10% to the total 
sum energy). We find that the total energy similarly de- 
pends on the spin to the result for a normal neutron star. 
The amount of the emitted energy increases as a spin is 
large negative (anti-parallel) , on the other hand it will be 
smaller for a large positive (parallel) spin. This result is 
similar to that for a normal neutron star, but the reason 
is not simple in this case. In the case of a normal neu- 
tron star, we could easily explain its spin dependence by 
considering a turning point. Because the location of the 
turning point monotonically changes with the spin (Fig. 
H) , then the gravitational waves would also be more radi- 
ated near the turning point as the spin gets large negative 
(anti-parallel). However for the case of an ultracompact 
star, from Fig. g we find that the location of the turning 
point has a minimum. Hence the change of the turning 
point does not explain the monotonic spin dependence of 



the total energy. We should point out that since the par- 
ticle passes through very strong gravitational field, the 
coupling between an orbital angular momentum L and a 
spin 5* is not simple, b ut m ay depend on the higher order 
terms of M/r in Eq. (3.1). 

For an ultracompact star, the ratio of the total angu- 
lar momentum to the total energy of gravitational waves 

As we mentioned be- 
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{Jgw/Eq-w) is shown in Fig. 
fore, the characteristic frequency of the system is impor- 
tant to think JcwZ-^'GW- However, for an ultracompact 
star many quasi-normal modes are excited by a scattered 
particle. Therefore, we must roughly estimate (-Eqw); 
and (Jgw)/ by using two frequencies: One is the or- 
bital frequency, and the other is the quasi-normal mode 
frequencies. Using these two characteristic frequencies, 
(-Eqw)/ and (Jgw)/ can be described formally as 



(£^Gw); 



(axial) 



(s(™"*)),+;^(4^™') 



(Jgw) 



(axial) 
I 



(j(= 



^)/+E(-^^™^)'' 



(3.6) 

(3.7) 



where (ii^^'^""*)); and (J^'^""*^); are the total energy and 
angular momentum from the continuous spectrum, and 
(Sj ); and (Jj ™ )j are ones due to the z-th quasi- 
normal mode. We expect that each contribution has the 
same relation (|3.5|), i.e., 




(axial) 



(axial) 



1 



^(cont) ' 



(qnn 



(3.8) 
(3.9) 



(qnm) 



Since ujI is independent of S, we believe that the 

spin dependence in Fig. |l^ is only from the continuous 
spectrum, whose characteristic frequency contains two 
parts: One is the orbital effect which corresponds to the 
angular velocity at the turning point and the other is the 
strong gravitational field effect which is excited to the 
motion of the particle. 



IV. CONCLUDING REMARKS 

We have studied axial modes of gravitational waves 
from a spinning particle scattered by a uniform density 
star. Particularly, we have focused on a spin effect to the 
energy spectrum, the waveforms, the total energy and 
angular momentum of gravitational waves, and analyzed 
a dependence of the emitted gravitational waves on a 
particle spin. 

In the case of a normal neutron star model {R ^ 3M) 
to which an observable neutron star belongs, we find that 
the total energy oi I — 2 gravitational waves gets larger 
as the spin increases in the anti-parallel direction to the 



orbital angular momentum. For higher multipole modes, 
however, it is found that the total energy does not de- 
pend monotonically on the spin, but has a maximum. 
The energy spectrum have only one broad peak induced 
by the particle encounter. This broad peak is affected 
by a spin of a particle, because the trajectory is deter- 
mined not only by the energy E and the orbital angular 
momentum L, but also by the spin 5*. In fact, we find 
that the broad peak of the energy spectrum shifts to the 
higher frequency region, as the spin gets large negative 
(anti-parallel). This is because the angular velocity of 
a particle at the turning point increases monotonically 
as the spin gets large negative (anti-parallel). However, 
it is known that this tendency of the broad peak shift is 
also found by a change of the orbital angular momentum. 
Therefore for a normal neutron star model, it is very dif- 
ficult to distinguish a spin effect from the orbital angular 
momentum effect in the energy spectrum of axial modes. 

For an ultracompact star model {R ^ 3Af ), in addition 
to such an orbital contribution, we find sharp peaks in 
the energy spectrum, which completely correspond to the 
quasi-normal modes excited by the scattered particle. If 
a particle passes through very close to an ultracompact 
star (rmin ^ 3M), the trapped quasi- normal modes give 
a prominent contribution to the total energy of the emit- 
ted gravitational waves. For example, in an ultracompact 
star with R — 2.26M, first eight trapped quasi-normal 
modes dominate. In particular in the case of Fig. Q the 
contribution from fifth to eighth peaks gives about 45 ^ 
70 % of the total energy. But if the turning point is far 
from a ultracompact star, the contribution from trapped 
quasi-normal modes is less than 10 % (Fig. ^. Note that 
in the case of a particle passing through very close to an 
ultracompact star, the frequency of a broad peak shifts 
monotonically as a spin increases, however the turning 
point and its angular frequency do not change monoton- 
ically and have the minimum value at some value of the 
spin. Therefore the interpretation of the cause of a broad 
peak is not simple. This unusual dependence of the turn- 
ing point on the spin might be caused by a complicated 
coupling between an orbital angular momentum L and 
a spin S, and this coupling may appear when a particle 
passes through in the strong gravitational field near the 
maximum of Regge- Wheeler potential. 

With respect to the ratio of the total angular mo- 
mentum Jgvv to the total energy -Eqw of the emitted 
gravitational waves, we could understand that its depen- 
dence on the spin is determined by the characteristic fre- 
quency, i.e., the angular velocity at the turning point 
or the frequency of quasi-normal modes. If the particle 
passes through the strong gravitational field, the ratio is 
also depend on the excitation of the strong gravitational 
field. Since the quasi-normal modes do not depend on 
the properties of a particle, the dependence of the ra- 
tio Jgw / Eqw on the spin is due to the change of the 
trajectory. 

From our results, we may conclude that the spin ef- 



fects on the emitted gravitational waves are seen only as 
the orbital change of a spinning particle due to a spin- 
orbit interaction in Eq. (2.1) and Eq. (p^). However, 
we expect another spin effect, which is a direct change 
of the source term of gravitational waves via a sp in con- 
tribution on the energy-momentum tensor ( 2.14 ). This 
may affect the emitted gravitational waves and its change 
could extract a characteristic property of a spin from the 
observed gravitational waves. However, in the case of 
a spinning test particle, our results show that the spin 
effect is mostly found as the orbital change by the equa- 
tions of motion. In fact, we expand the source term with 
the value of the spin s for the first order in order to exam- 
ine the particle spin contribution to the spectrum from 
the energy-momentum tensor. We compute the spec- 
trum in the case of the source term with no spin, with 
first order spin, and with full order spin. There are no 
remarkable difference among these three results, and all 
of them has a global peak near the turning point for the 
normal neutron star. Therefore, we conclude that the 
spin effect via the spin part of the energy-momentum 
tensor is much less than that by the orbital change. In a 
realistic binary neutron star, however, a spin part of the 
energy-momentum tensor plays definitely an important 
role in the gravitational wave, because it has a stellar 
structure and it is essentially different from that of the 
point particle, and the emitted gravitational waves may 
depend on such a term as well as the orbital change. 

In a realistic situation, the polar modes are also impor- 
tant, because those contain the fiuid oscillation modes 
and then /, p modes have some information about high 
density fluid matter of a neutron star. Then the effect 
of a spin on those modes may reveal more information 
about neutron star matter. Although we could analyze 
those modes in the present model, that is a perturba- 
tion approach with a spinning test particle, we expect 
that the results will be negative. Because the different 
interaction between a particle spin and polar perturba- 
tions from axial perturbation will appear mainly in the 
energy-momentum tensor. However, as mentioned above, 
since the contribution from the energy-momentum ten- 
sor would be very small in the perturbation approach, 
we believe that the result for the polar modes may be 
similar to the axial case, i.e., the main factor of the spin 
effect is caused by the orbital change and the difference 
of /, p modes due to the particle spin appears only in 
an enhancement of its strength. On the other hand, if a 
star is rotating, we expect a coupling between a stellar 
rotation and a spin of a particle, and then it is important 
to analyze it, particularly because a rotating star shows 
the r mode instability. 
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APPENDIX A: SOURCE TERM OF A SPINNING 
PARTICLE FOR THE AXIAL MODE 

For a spinning particle, the difference in the analysis 
from our previous work (non-spinning particle) [ |lO| ap- 

,(cxt) 



pears only in the source term 5;,,^^ in Eq. ( 2.2(: ). Since 
the detail of our method is shown in Ref |10||, we only 
present the part of the calculation concerned with the 
source term of a spinning particle. 

The amplitude of the gravitational wave is given by 
Eq. (A24) in Ref |1^ and then we must calculate the 
following integration. 



dr*xt^:^^'\r*)St:^{r*). (Al) 



When a spinning particle is scattered on the equatorial 
plane by a relativistic star or a black hole, with its tra- 
jectory as z^{t) — (t{T),f{T), 9 — tt/2, ^(t)), we find Eq. 
(Al) as follows: 
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Pirn is the associated Legendre function, and QmL , QmL , 



i(o) nil) 



DinL, D\nL and D)nL are given as 
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Ks is written as 
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where n ^ l{l + l)/2 — 1, and Cim is a normalization 
constant of a spherical harmonics Yim , written as 

Qm = (-D'^+i-i^/^/^ ^^ilTl^ii; . (A3) 

y in [l + \m\y. 



which is defined in the relation between the 4- velocity v^ 
and the normalized momentum u'^ as 
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FIG. 1. Timelike condition for the 4- velocity and the scattered orbit condition for a spinning particle are shown in the case 
of (a) S = O.bfiM and (b) S — O.SfiM. The dotted line denotes a null condition, i.e., d''u^ = 0, and the timelike condition 
is satisfied in the region above the dotted line. The solid line corresponds to the maximum point of the effective potential 
y (particle) rpj^^ shaded region presents a forbidden parameter for a scattered orbit. 
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FIG. 2. Location (a) and angular velocity (b) at the turning point of a spinning particle in the case of -R = 5.0M , E = l.Ol/j., 
and L = A.b^iM. 
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FIG. 3. Location (a) and angular velocity (b) at the turning point of a spinning particle in the case oi R = 2.26M, 
E = 2.30^, and L = 12.0^^M. 
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FIG. 4. Energy spectrum {1 — 2 mode) of gravitational waves from a spinning particle with E — 1.01/i and L — 4.5/iM 
scattered by a uniform density star with R — 5.0M. In the parallel case (a), solid, dotted and dashed lines show the case of 
S = 0, 0.2/iM, and 0.8/iM, respectively, while in the anti-parallel case (b), solid, dotted and dashed lines show S = 0, — 0.2/iM, 
and —0.8/iM. Vertical lines represent the angular velocity of a particle at the turning point. 
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FIG. 5. Total energy of gravitational waves for each multipole mode I in the same situation as Fig. |^ Solid, dotted, dashed, 
dash-dotted and dash-two dotted lines show the case of i = 2, 3, 4, 5, 6, respectively. 
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FIG. 6. Ratio of the total angular momentum to the total energy radiated as gravitational waves for each multipole mode 
I in the same situation as Fig. 0. Figure (b) is the ratio (a) multiplied by the angular velocity at the turning point for each 
spin value, in order to estimate Eq. (3.5). Solid, dotted, dashed, dash-dotted and dash-two dotted lines show the case of 
/ = 2, 3, 4, 5, 6, respectively. 
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FIG. 7. Waveform {1 = 2 mode) of gravitational waves in the direction oi d = 7r/2, = in the same situation as Fig. ]A. 
In the parallel case (a), solid and dotted lines show the case of S = and 0.8/iAf, respectively, while in the anti-parallel case 
(b), solid and dotted lines show 5* = and — 0.8/iAf. 
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FIG. 8. Energy spectrum [1 = 2 mode) of gravitational waves from a spinning particle with E = 1.01/^ and L = 4.5fiM 
scattered by a uniform density star with R = 2.26M. In the parallel case (a), solid, dotted and dashed lines show the case 
of the spin S = 0, 0.2/^A4", and O.SfiM, respectively, while in the anti-parallel case (b), solid, dotted and dashed lines S = 0, 
— 0.2/^M, and —O.SfiM. Vertical lines represent the angular velocity of a particle at the turning point. The small figure in (a) 
and (b) is the enlargement of 1st quasi-normal mode, in order to show the difference of its strength. 
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FIG. 9. Energy spectrum {1 = 2 mode) of gravitational waves from a spinning particle with E = 2.38/i and L = 12.0fj,M 
scattered by the same star as Fig. H (i? = 2.26A4). In order to focus on the detail of peaks in (a), we show the enlarged figure 
(b) for the frequency range of 0.0 < Mu < 0.5. Solid, dotted and dashed lines show the case of the spin S = 0, 0.2/iAf , and 
0.8^M, respectively. Vertical lines represent the angular velocity of a particle at the turning point. 
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FIG. 10. Location (a) and angular velocity (b) at the turning point of a spinning particle in the case of 7? = 2.26M, 
E = 2.38^, and L = 12.0/xM. 
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FIG. 11. Waveform {1 = 2 mode) of gravitational waves in the direction oi 9 = it/2, <j) = in the same situation as Fig 
Solid and dotted lines show the case of the spin 5 = and 0.8/iM, respectively. 







FIG. 12. Total energy of gravitational waves from a spinning particle with E = 2.30/^ and L = 12.0/iM scattered by a 
uniform density star with R = 2.26M for each multipole mode I. Solid and dotted lines show the case of / = 2, 3, respectively. 
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FIG. 13. Ratio of the total angular momentum to the total energy radiated 1 — 2 gravitational waves in the same situation 
as Fig. ^. 
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